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ABSTRACT 

The  stability  of  a  toroidal  gas  discharge  is 
examined  by  using  an  energy  principle.   In  the 
equilibriiam  configuration  the  plasma  is  assumed 
to  be  perfectly  conducting  and  to  have  constant 
pressure  and  density.   It  is  found  that  the  pinch 
can  be  stabilized  under  certain  very  restricted 
conditions.   However,  one  of  the  conditions  Is  that 
the  aspect  ratio  be  small  so  that  a  long  thin  torus 
will  always  be  unstable  in  the  limit. 
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THE  HYDROMAGNETIC  STABILITY  OP  A 
TOROIDAL  GAS  DISCHARGE 

1.   Introduction. 

In  recent  years  many  experiments  on  the  pinch  effect 
have  been  carried  out  using  both  cylindrical  and  toroidal 
geometries.   Until  now  cylindrical  gas  discharges 
have  been  treated  theoretically  and  it  has  been  found  [IJ 
that  the  pinch  can  be  stabilized  under  certain  conditions, 
It  is  the  purpose  of  this  investigation  to  study  the 
stability  of  a  toroidal  gas  discharge  and  to  get 
some  indications  of  the  influence  of  the  curvature  on 
the  stability.   This  Influence  should  be  not  negligible, 
if,  as  was  the  case  in  some  experiments,  the  so-called 
aspect-ratio  —  the  ratio  of  the  small  diameter  of  the 
torus  to  the  large  diameter  —  is  not  small  compared 
to  one. 

Ti-jo   different  methods  may  be  applied  for  the 
Investigation  of  the  stability  of  a  hydromagnetic 
equilibrium  configuration.   The  first  method  makes  use 
of  the  normal  mode  analysis  and  with  this  method  the 
stability  of  cylindrical  gas  discharges  have  been 
treated  [IJ.   But  for  more  complicated  equilibriixm 
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conf igiirstlons  it  tiu-'ns   out  that  the  application  of  an 
energy  principle  is  more  useful.   This  method  is  used 
in  the  following  treatment «   The  main  reason  for  this  is 
that  the  mode  analysis  leads  to  a  differential  equation 
where  the  boundary  conditions  are  quite  com.plicated  in 
the  case  of  toroidal  geometry.   This  complication  can  be 
avoided  using  the  energy  principle. 

2*   The  Equilibrium  Conf iguration» 

We  assume  that  the  plasma,  regarded  as  an  ideal, 
compressible,  perfectly  conducting  fluid,  has  the  shape 
of  a  torus  and  is  surrounded  by  a  conductor,  as  indicated 
in  Figure  1,   In  the  equilibrium  configuration,  the 
pressures  and  density  inside  the  plasma  are  taken  to  be 
constant.;  it  is  assuraed  that  there  are  no  electric 
currents  flowing  except  on  the  surface  of  the  plasma; 
and  the  m.agnetic  field  is  taken  to  be  of  the  form 

B  =   (0,b/r_,0)   inside  plasma  (2.1a) 

-^ 

B      =      (B   ,b/r,B    )      outside   plasma  (2.1b) 

in  which  the   components   have   been  listed   in  the    order, 
r-component,      ^-component,       z-com-ponent ,      and    in  which 
b      and     b      are    constants.      The    coordinate    system 
(r,^,  z)      is    cylindrical,    as    indicated   in  the   fig^areo 

In  order  that  this  be   an  equilibrium,   configuration. 
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the  meridional  external  field,  v;hose  components  are 

A.  -A 

B   and  B  ,   must  (jointly  with  the  given  azimuthal 
fields)  provide  the  required  confining  pressures  on  the 
plasma.   The  field  must  be  purely  tangential.   The 
pressiore  balance  gives  the  relation 

I  (B^  -  B^)   =  p  (2.2) 

where   p   is  the  constant  equilibrium  pressure  of  the 
plasma.   This  external  field,  if  extended  in  the  stirround- 
ing  space,  would  have  various  singularities,  for  example, 
on  the  axis.   It  must  therefore  be  supposed  that  somewhere 
outside  the  apparatus,  but  not  shown  in  the  drawing,  there 
are  suitable  current  coils  to  produce  this  field.   The 
conducting  wall  shown  is  supposed  to  have  no  effect  on 
the  equilibrium  field.   The  cross-sections  of  the  plasma 
surface  and  of  the  conducting  wall  are  assumed  to  be 
circles  lying  on  coordinate  surfaces  \    ~       \.      ^J^d. 
^      =  \^      in  a  toroidal  coordinate  system   (^,^,9) 
and  therefore  not  quite  concentric.   With  this 
assumption,  the  conductor  is  not  placed  on  a  line  of 
force  of  the  vacuum  magnetic  field  as  the  equilibrium 
condition  requires.   But  we  may  assume  that  the  equili- 
brium field  has  enough  time  to  penetrate  the  conductor 
and  that  the  conductor  acts  only  on  the  perturbed  field. 
This  Inconsistency  and  the  perhaps  improbable  nature  of 
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the  current  coils  that  would  be  required  to  produce  the 
equilibrium  field  specified  [L|.l,  are  excused  on  the 
ground  that  the  finer  details  of  the  equilibrium  field 
probably  have  little  effect  on  the  hydromagnetic 
stability  and  would  not  be  accurately  known  in  a  given 
experimental  set-up  anyway. 

The  toroidal  coordinates  are  related  to  the 
cylindrical  coordinates  by  the  transformation 

_     as  inh  V?  r  o  o  „  ^ 

cosh  >i  -  cos  9 

0=0  (2,3b) 


a  sin  0 /^  ^^\ 


^     cosh  h  -    cos  9 


where      (+aj,0)      are  the   poles   of  the    coordinate    system 
in  the   x,y~plane.      Prom  this  we   get   the   radius    of   the 
'cross-section  of   the   plasma 

Bo  =      TTStTTT  ''•'^' 

■o 

with  its  center 

cosh  ^ 

r  =   a  — r— ; (2.5) 

o      Sinn  M 

and  the   radius   of   the    cross-section  of  the    conductor 


R,    =        .    ,  ^  ,,         .  (2,6) 

1  smh     >{-. 
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¥e  define  the  asDect  ratio 


^1 


^  =  _i  =  cosh  ^^  .  (2.7) 

o 

The  equilibrium  condition  (2o2)  is  given  in  this 

coordinate  system  by 

-I   r  ^o      (cosh   »?   -COS  e)^    A-3    9   / 

l)(g2+_^ (b2-b2)    =  p.  (2,8) 

-^^  I  ^     a^  sinh^  Yj^  J 

3o   The  Energy  Principle. 

As  mientioned  above  we  want  to  use  an  energy 
principle  [2 J  for  investigating  the  given  equilibrium 
configTxrs.tion.   This  energy  principle  says  that  the 
necessary  and  sufficient  condition  for  stability  is 
that  the  change  in  energy   5M   for  a  given  displacement 
t,     must  be  positive  for  all  £,    o 

The  change  in  energy  has  been  derived  by  Bernstein 
et.  al.  [2]  and  is  given  by 


5W     =      5W 


(i)    +   g.^(e)    +   g^(s)    ^  (2^-^^ 


/    .     A 

Here      5¥  is    the    contribution  from   the   plasma    (the 

internal  part) 


5w(i)      =     i 


r    -^     •?>  ^ 


[Q  °Q  +  YP    (div  C)^J    dV  (3o2) 


assuming  that  the  pressure  is  a  constant  and  that  no 
electrical  currents  are  flowing  in  the  plasma.   The 


volume    Integral  has  to  be   extended  over   the    initial 
volurae   of  the   plasma,      Q,     is    given  by 

Q     =      curl    (^  X  B)  (3.3) 

and     Y      is    the   ratio   of   the    specific   heats.      The    contri- 
bution  of   the    vacuum      5W  ®         (the    external   part)    is 


given  by 


5W^^^      =     I     /(grad  1)^   dir  (3,1+) 


where     J     is   the    scalar  potential   of   the   perturbed  vacuvim 
field  and   the    integration  must  be   extended   over   the   vacuum 
which   is    surrounded  by   the   plasma   and  the    conductor. '    The 
potential      ^     must   fulfill   the    two    boundar*^    conditions: 


^     =      curl    (C  X  B) 


(3.5) 
n 


on  the   plasma-vacuum    interface,    where    the    subscript   refers 
to  the   normal    component   and 

on  the   containing  wallc      (^J/t^n     is   the   normal   derivative 
of     ^,      These    two   boundary    conditions    do   not    determine 
the   potential   uniquely   in   some    cases,    since      2!     need  not 
be   a   single-valued  function  necessarily    [3]'      This  point 
will      be   discussed  later.      Finally,    the    contribution 
from  the   plasma-vacuum  interface   is   given  by 
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5VJ 


^^      ■■=     \    J  (n°4)^  (no-(grad(p+  |!B!2)>)da-     Oo?) 


where  <l .  ^  ^^     denotes  the  increment  of  a  quantity  across 
the  boundary  in  the  normal  direction  no   The  integration 
has  to  be  extended  over  the  plasma- vacuum  interface. 

Using  the  energy  principle  the  equation  of  stability 
is  reduced  to  an  examination  of  the  sign  of   5¥  for 
arbitrary  displacements   ^<,   To  do  this  one  must  minimize 
5V/a   For  this  procedure  we  follow  a  method  suggested  by 
Bernstein  at.  alo  [2 J   First   5W^^''   and  bW^'      are 
minimized  separately  by  prescribing  n  "  E,      on  the 
plasma-vacuum  boundary. 

Minimizing  6W     leads  to  the  Euler- Lagrange 
equation  for   C 

Y  P  grad  (div  C)  -  3  x  curl  Q  =  0  ,       (3«8) 

When     C      is    chosen  to   satisfy    (3o8),    the    integral 

SW  in    (3" 2)    can  be  reduced  by  partial    integrations 

to   a   surface    integral, 

"^^^        Wn     =|J(n°€)    ^rp   div  C-B«Q]    d  a-    .  (3<.9) 

(e )  -r 

The  external  part   5W •  '   is  minimized  by  choosing  ^ 

to  satisfy  the  Laplace  equations 

Ai   "    °  '  (3oio) 

When     J     is    so    chosen^,   we   get 
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[5W 


(e) 


W  "  I  /^^  ''■'  S^^^  ^^  ^  °^*      ^^'-^^^ 


The  unit  normal   n  points  into  the  plasma  here,  and 
out  of  the  plasma  in  (3. 9). 

In  this  way   5W   is  now  reduced  to  a  surface 
integral  over  the  plasma-vacuimi  interface.   The  process 
of  minimization  is  completed  by  minimizing   5¥  given  by 
(3.1)  and  (3-7),  (3.9)  and  (3.11)  with  respect  to 
(n  o^)  . 

l|o   The  Stability  of  a  Cylinder. 

Before  we  calculate   5W   for  a  torus  we  will  show 
how   5W   can  be  obtained  for  a  cylinder.   In  a  later 
section  we  will  compaj^e  the  results  of  the  cylinder  with 
those  of  the  torus. 

In  a  cylindrical  coordinate  system   (r,j2^,z)   the 
.  equilibri'um  field  is  given  by 

B  =   (0,0, b  )       insiae  plasma         (I|..  la) 

B  =   (0,Bg(RQ/r,b  )   outside  plasma        (Li-.lb) 

and  the  equilibriiom  condition  is  given  by  the  pressure 
balance 

i    (^|-^^°c  -  ^c^  =     P  '  ^^-2) 

(In  bending   the   cylinder   into   a   torus,    the    z-component 
of   the  magnetic   field  becomes    the   ^-component    and   the 
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jZ5- component  becomes   the  meridional   component.      See 
eq.,     (2»l)    and    ([(-.l).-)      £,      can  be   Foiirler   analyzed   in 
the    form 


^     =     ^(P)   3i(a^+kz)      ^  (^^3) 


Now   the   problem   splits    into    two  branches    according 
as      k  7^   0      or      k  =   0    . 

(a)      k  7^  0    : 

Prom  (3«8)  it  follows  from  scalar  multiplication 
by  B   that 


Introducing 


(3,8)  gives 


and 


div  4   =   0   .  (1^.1+) 


^  =   -  (B  .Q)  /  bf  (l+o5) 


A^  =    0  (^-6) 


•*■•>■  •*■ 

(n  "  4)   =  — 5-  (n  o  grad  "J'  )  .  (^.7) 

k^ 

Using  (i^.,5)  and  (I4..  7)  we  get  for  (3.9)  s 

^^^^^^W  "  l~l  f$  (n.grad^  )  d  cr.        (i|.c8) 

For  calculating   5W^^^  (eq.  (3.7)),  the  differences  of  the 
total  pressure  gradient  across  the  surface  is  needed* 
Using  (1^., la)  and  (l^.lb)  we  get 
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(n  »  /grad(p  +  |  IbI^)^)   =   -  b|  /R^  ,  ()|.9) 

The  solutions  of  the  differential  equations  (Li.<,6)  and 
(3«10)  for  '^    and  J  respectively  are.  In  cylindrical 
coordinates: 

^    =     1    (kr)  e^^'^^"'^^^  (1^.10) 

-•-a 

and 

2  =   1       1       (k  R     ^     +   ah    )   i'(kR    )   »M(r)ei^'^^-*-^^'^    (l^oll) 
^pjj^2ocj3ao  ^ 

o 

where 

l'(kR-,)K   (kr)-K' (kR,  )I    (kr) 

M(r)      =  -,^       ^    ,^ -? ^--^ .  (ij..l2) 

I^(kRjK^(kR^)-lJkR^)KjkRj 

I    and  K   are  the  Bessel  and  Hankel  functions  of  the 
a        a 

»        t 

first  kind  and  of  the  a-order,   I„   and  K^   are  the 

a        a 

derivatives  of   I    and  K  ,   The  constants  of  integra- 
tlon   in  (Lj.,  11)  are  chosen  in  such  a  way  that  the  two 
boundary  conditions  i3'S)    and  (3«6)  are  satisfied.   In 
cylindrical  coordinates  they  have  the  form: 

~i-  =  -ikb  C   -  J^  4r  U  (2p+b^-b^)  (for  r  =  R    (11,13) 


and 


-T—  =  0        for  r  =   R,  .  (I)_.ll|) 
Finally,  from  (l+olO)  and  (liell)  we  find  for   5W 
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5W  =  -^(i'{kR    ))2    f(kR    )V   V^  + 
o  a        o 

where      L      is    the    length   of    the    cylinder.      The    first    tvjo 
terms    are    the    contributions   from   the   plasma   and  the 
vacuiira  respectively   and  they    are    always   positive,   while 
the    last    term   is    the    contribution   of    the    surface    and    it 
is   always   negativec      This    dispersion  relation  has  been 
given  for   instance  by   Taylor   [ij   for   given  values   of 
k      and      a.       The    sign   of      5W      depends    on   three   parameters 
which  might   be      b^  /B^s      b  ,  /B^      and     R^  /R^    . 

(b)      k   =  0 

Prom  (3o8)  it  follows  that  if   a  7^  0 

div  C  =   (B  «  Q)   =  0  ,         (l^cl6) 

Therefore  according  to  (3° 9) 

5¥^^^   =   0  .  ()+.17) 

The  solution  of  (3-10)  for  2^  which  satisfies  the 
boundary  conditions  (i|..13)  and  (liolL^.)  is  given  by 

o  1  o 
where   C  (R  )   is  an  arbitrary  constant  which  might  be 
normalized  to  unity <,   Using  ([|.ol8)  we  find  for   5W,   if 
k  =  0   and   a  7^  0 
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We  see   that   5W  >  0   for  all   o-  >  1   as  long  as 

R  /R,  7^  Oc   If   R  /Rt   approaches  zero   5W  >  0   for  all 
0   1  c  1   -  ^ 

a  >  1   and   5¥  =  0   for   a  =  lo 

The  case   a  =  0   has  to  be  treated  separately. 
From  (3e8)  it  follows  that 

4p  =  ar  (i4-o20) 

where      a     is   an  arbitrary   constant.      Prom  the  boiondary 
condition    (14.=,  11)    for     J^     we   get-^ 

■T-     =      0  for     r   =  R     .  (i;,21) 

Therefore 

5W^®^      =      0        .  (1|.22) 

Using    (l+clS)   we    find  for      5W     if     k  =   a  =  0: 

5W      =      TiL  a^   R^    r2p(Y-l)+b^+b^V     «  (l;o23) 

01  c      cj 

Here  SV/  is  always  greater  than  zero  since  \  >_  1, 
Therefore  the  cylinder  is  com.pletely  stable  against 
radial  oscillations  as  one  would  expect. 

Prom  (i^<,l9)  and  {[j.o23)  we  can  conclude  that  the 
cylinder  is  completely  stable  against  all  perturbations 

— /  In  this  case   J   is  not  a  single -valued  function.   But 
J  does  not  contribute  to   5¥   according  to  (li.,20)  and 
therefore  this  is  not  of  ijnportance  for  the  cylindero 
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which  are  independent  of   z   with  the  exception  of  the 
case   a  =  1   if  there  is  no  conductor  present.   In  this 
case  the  equilibriinn  is  neutral  for  a  perturbation 
a  =  1, 

5o   Calculation  of   6W  for  a  Torus. 

(a)  The  solution  of  Laplace  equation  in  toroidal 
coordinates. 

In  toroidal  coordinates  E,      can  be  Fourier  analyzed 
in  the  form 

t    =   €  (>^,6)  e^^^^-"^^^  ,  (5ol) 

The  problem  now  splits  again  into  two  branches  according 
as  m  7^  0  or  m  =  0,   since  m   corresponds  to  k  in  the 
cylindrical  case.   The  case  m  =  0  will  not  be  discussed 
in  this  paper.   A  close  investigation  has  shown  that  the 
mathematical  formulation  is  more  complicated  than  for 
the  case  m  7^  0   and  in  addition  2!j   the  scalar 
potential  of  the  vacuum  field,  is  not  always  single-valued 
as  it  is  when  m  7^  0.   However,  from  the  nature  of  the 
perturbations  and  from  the  results  derived  above  for  the 
cylinder  when   k  =  0  we  expect  that  with  an  outer 
conductor  the  torus  will  be  completely  stable  against 
perturbations  with  m  =  0, 
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In  the  case  m  ^   0 ,    the  equations  (I4-.I4-)  to  (I4.08) 

hold   also   for   the   toroidal    freorcetry    if  we   replace     b 

c 

and      k     by     b      and     in     respe  ctively  o      Pur  the  more   using 
(2.1a)    and    (2olb)   we    get 

(n  o(grad    (p   +  |    \b\^)})      =     -  A    [2p  +  ^   (b^-b^)]  .   (5.2) 

o  r 

In  this  way    (L|.o8),    (3<,ll)    and    (3o7)   become 


6W(1'      =      -K^ 


m 


(^Tt 


o 


r   \!/  4^  d© 


(5o3) 


5W^^^^      =      71 


n2TZ 


r  i  -y-  de 
o  ^\ 


(5.1;) 


and 


5W 


(s) 


a 


.2% 


Jo 


(C^)^[2pr2-f  ^   (b2»b2)]    d9      (5.5) 


where      5        is   given  according   to    (l+^y),    by 


a   sinh     v^ 


HI 


(cosh  V|   -   cos   e)      '^^ 


(5c6) 


The  two  boundary  conditions  {3»S)    and  (3.6)  become 


A 

r    sinh  Vi      A        r 


^Tr   €^[2p+^    (b242)]l/2J    (5„7) 


for     y\  =     ^        and 


=    0 


(5.8) 


for    1    =     1^. 
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The  functions  ^  y  J  appearing  in  (5«3)s  (5ol4.), 
and  in  (5c  6),  (So?)  are,,  as  we  have  seen,  solutions  of 
the  Laplace  equation  and  therefore  may  be  made  up  of 
sums  of  terms  of  the  type 

■/cosh  y|^"COS  e  JA^  q"^  3_  (cosh  )^)+B^P"'  ^  (cosh  ^ )  p^"'^'^®^ 

where   P^  ,  and  Q   -  are  the  associated  Le^eendre 
1       „  1 

functlonso   As  ^  is  regular  at   W  =  oo  ,   the  axis  of 
the  discharge^  it  will  consist  of  Q-functions  alone; 
we  shall  choose 


$j^  =  /cosh  V|  -  cos  0'  Q,"^  ^  (cosh  r\)   e  (5<.9) 

""2 


where  here  and  in  what  follows  the  argument  of  the 
Legendre  functions  as  well  as  the  ^-variable  have  been 
suppressed  for  simplicity  in  writingo   Then  the  corres- 
ponding exterior  function  will  be 

■^n     ^      (        ^^    s  3„i   s  ^„1 

The  constants  ^a  ,  '^''b    are  to  be  determined  by  the 

s*    s 

boiondary  conditions  (5°  7)  and  (3o8)o   Clearly  the 
functions  (5«9)  form  a  complete  set  of  allowable 
functlonSj,  but,  unfortunately,  these  functions  are  not 
orthogonal  with  respect  to  the  mode  number  n   (clearly 
they  are  orthogonal  wi1-"      -ot  to  m   and  various 
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m-modes  may  be  separated).   This  raeans  that  one  must 
investigate  not  only  all  n-m.odes,  but  also  all  linear 
combinations c   Thus  one  must  form  the  integrals 


where 


""   nf  2m2    Jo  ]^^  W       ^^  ^/ 


f  f^kc),  2|-i:^^i^j^^ 


5W^|^      =     5      1         r<)T„  c)«    ?^'  +  ?!'().■$.(  d@  (5.12) 


1271 


2        r    (b2-^2) 


_<fr 


the   star   denoting   complex   conjugate «      The   necessary   and 
sufficient    condition  for    stability   is    that    the      matrix 
(6W    />)      be   positive    definite    (or    zero)    for    all     m<, 

(b)    The   Plasma-Vacuum  Boundary,, 

In  this   section  vie  work  out    explicitly  the  boundary 
condition  of   equation    (5»7)o      Substituting    (5^6)    into 
(5o7)    gives 

^,ii|,asinh^l  „osh„-cos  9  );^f ^^-^    f ](S.13) 

A 

where   B„   is  the  vacuum  value  and,  according  to  the 
y 
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pressvtre   balance   equation   (2o8),    is  given  by 


"^a   sinhT   r-    cosh  y(_-2   cos   9)"   +  M"^   sinh^v^^ 


2   "^2 
It   is   assumed  that      b   -b     >   0,      which  is    the   usual    case. 

It    can  be    shown   that    all    our   results    are    correct   also 

when     b      >  b    ,,      provided   only    that    one    Interchang-e      b 

/\ 
and     b      so    that   the    square   roots   in    iSollO    and    (5.21) 

2 
remain  real   and     M        positive..      When    (5o9)    and    (5=10) 

are   employed  we    obtain  the  results 

^^n     ._     sinhn  j       «    i(n-l)e,/       ,,    f        ,s    inG 

^.^     _     — __— __ 4__ — _  ^  »q  e    ^  +(q  +2q    coshh.ie 

MT  2  Vcoshvi-^cos   Q    I    ^  ^     ^  ^ 


2  y  cosh  Yi  -cos   © 

^»      i(n+l)e 
•=   q     e 


(5oi5) 


and 

^^     ._.  sinh  h  ^  n.        isG  ic  i^    ) 

I-  2vcoshv|-cos    ©      s==oo 

def  „  , 

s  s~l^s-l        s-^s-1  s^s        s^s 

(5.16b) 

+  2   cosh      (^a^q'+"b^p' )  -  (^a^,.q'    ,+^    ,,?  '     )  . 
s^s        s-^s  s+l^s+1       s+l-^s+1 

def     m  def  m 

Here  we   have    set,    for    short      q        =Q       t,P=P       tI 

'  ^n     ^   1  '  -^n      1  ^ 

the  prime  denotes  differentiation  with  respect  to  the 
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argument,      cosh   /|  »      The   Index      "n        in  the   upper 
left-hand   corner   indicates   that   we   are    fitting  to   the 
n  plasma  mode   as   given  by    (5''9.)« 

The  procedure  now    is   to   set    (5ol5)    into    (5ol3) 
and   to    expand  the    right-hand   side    as    a   Poiirier    serieso 
This    leads    to    the   Fourier    coefficients 

^C      =  ^  ^G     +  .a-=^ii%JL    f       (cosh  V,"COs   e)3/2 

R  1  ^5.17) 

e  .^—  n^    ^ie(  ^~ise   .^ 


(coshh-cos   0)   ■'        J 


where        G        are   the   coefficients   of   the   3-term 
Fourier   series    in  the  brackets   of    (5'-l5)»      A  partial 
integration  may  be   performed   on   the   right   hand   side 
of   this   expression  and  we    obtain 

r,        -      o   b   n_  s  Inhh   -/,  2  C2  v—  n^      ^        c  r-  i  o  i 

C        -      1  —     G^    -  ■ —      o  -  7b    "b     >        G.    I.„    (5olo; 

n  s  m        s  OTrK„2      '  ^      /     /s 


2Trm 


where      I^        are   the   integrals 


sin  9-is  (  cosh  >] -cos    Q)    f  p 

^  Ms   cosh>i-2   cos    e)"^ 


2 
(cosh  h  ~  cos  0) 


■^m2  sinh2^(^/2  ^i(l~s)e  ^Q 


(5.19) 


The  integrals  I^    can  be  evaluated  by  standard  methods 

and  we  obtain 

^/s  "  I  7iiMe"''^"^ho(3^5s)    f or  i  7^  s 

(5«20) 
I^^  =  27:1/(1  -  |)  » 
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Setting  this  result  into  (5ol8)  gives  the  following: 

^C   =  il[b"~y^^^  sinhV"(l-4M}J  [q  +2q'  cosh  yj] 

■=  -  -/b^-b^  sinh  >)  (Me~  ^)  q'  C 
m  C        n  I 

/  K  pi    ) 

+   sl^it/b-=b'^  v2n+5)  Me"^(q^+2q^^  cosh  )^) 


^ 

are   now  i 

"^a 

= 

^C 

and 

°A 

=: 

0 

G     =   i_£i§il4.  Vb^I^  (Me-^'^JL  (2n+5s}q  +6q'  sinh  m] 
—        i\vrv 

for   s  >  2  » 


Our  boundary  conditions  which  determine        a        and 


for     Y\^  =     ''l^  (5.22a) 

for      7  =     ^2  (5o22b) 

where   ^A   and  ''^C„   are  given  by  (5ol6)  and  (S«2l) 
respectively.   For  each  n  we  have  a  different  system 
of  equations..   In  practice  the  infinite  set  of  equations 
(So 22)  will  be  replaced  by  a  finite  set,  with  s   running 
from  -N  to   N-1,   say;  there  will  therefore  be  l\N 
equations  for  each  n. 

(c)  Calculation  of  the   5W«s. 

The  integrals  involved  in  the  explicit  calculation  of 
the   5W»s   given  by  equations  (5cl2)  can  be  evaluated  in 
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closed  formo      In  fact,    all    of   these    integrals    are    of 
the   form 

il    =        I       ^ ^ (5o23) 

^  J  ^      (cosh  >|^-    cos   Qr 

where      X,ij,      are    integers,      Ptjt  the  more,      X     may  be   taken 
positive   without    loss    of   generality,    as   the    substitution 
0   — ^    _   9     leaves    the   integral    invariant   and     I        =  I    . 
The    integrals   needed   are    for      [j,  =   1   to    5      inclusive    and 
are  i^rritten  out   explicity  below: 

1        sinh  ^  2        3.^2  ^    ^  sinh  ^ 

3        2lsinh>  ^^^     *!  sinhH 


I 


.  lg'^-^3)e-^t  ^   ISell!]  (5.21,) 

slnh        ii  sinh-^n. 

I J   =  ^^^    H        [(X+ii)(X+3.)(X+2)(?.+l)    + 

-^        Ixi    slnh-'  \ 

-■»]  -2  >t 

+  10U+k-){-k+3)U+2)      .  ^,         +  [^5(A+li)(A+3)        ^   g        + 

L  sinh    \ 

+   105(A+I;)     ^3    +  105    ^  I   ]  . 
slnh  >Y       sinh^  >|^ 

Setting  (So  9)  into  the  first  of  (5=12)  gives 

m 
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or 

2  2  r  -\^~'^]t^rj 

Here   5  /,  is  the  Kr one cker- symbols   Similarly,  setting 
(5»10)  and  (5ol6)  into  the  second  of  equations  (5"12) 
gives  the  result 

S^l^^'  =  71^  a  sinh  >;  HT"  (""a  q  *%   pj    C   e  '""  ^  '  °  (.5.26) 
n£  (o  -^  ^'    s^s    s^s   r 

Finally^  for  the  third  of  ev  -ol2),  we  get 

n^         1-^     V  a"   sinh^W  /  \   ^^  ^  ^ 


Um 


Yo' 


-21q„qK,^!l,f-<i*q,q^4-^i)J.         (5.27) 

It  has  been  postulated  that  bending  a  very  long 
thin  cylinder  into  a  torus  will  not  alter  the  stability 
properties o   That  this  is  not  so  is  evident  from  an 
examination  of  the  above  equations  in  the  limit  ^   — >   0* 
In  particular,  it  is  quite  easy  to  show  that  for  m  =  1 
the   (OjO)   and   (0<,+l)   terms  in  the  matrix  all  behave 
in  the  same  way  in  the  limit  and  thus  the  matrix  never 
becomes  diagonal  as  would  be  the  case  if  a  long  thin 
torus  had  the  same  properties  as  a  cylinder.   (This  result 
has  been  obtained  for   SW'-^'   and   5¥-^''c   For   6W '  ''   the 
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dependence  of  the  "constants"  ^a^   and  ^^   on  ^^ 
is  very  complicated  and  no  expansion  in  terms  h 
seems  possibleo   However,  the  niomerical  computations 

show  convincingly  that  the  contribution  to  the  matrix 

(e  ) 
elements    from      5W  behaves    in  the    same   manner   as    do 

the  other  contributionso )   Physically  this  means  that 

there  is  always  at  least  one  mode,  namely  m  =  1, 

which  feels  the  curvature o   The   m  =  1   mode  corresponds 

to  a  perturbation  having  just  one  wave  length  around  the 

torus  and  this  wave  length  is  always  of  the  same  order 

of  magnitude  as  the  curvat'ure. 

6.   Method  of  Computation  and  Results. 

Codes  have  been  x-;ritten  for  the  IBM-7OI4.  to  determine 
the  stability  of  a  plasma  for  both  the  cylinder  and  the 
torus.   For  the  cyliader,   SW   is  compu'^-'-d  directly 
from  eqo  (l4.ol5)<'   The  Bessel  functions  tb.at  are  needed 
are  evaluated  using  a  code  written  by  Goldstein  and  Kresge 
[5]  and  results  agree  with  those  publisher   /y  Tayler  [Ijc 
The  main  use  of  the  code  is  to  facilitate  comparisons 
between  the  cylinder  and  torus o 

For  the  torus,  the  situation  is  much  more  complicated. 
First  the  associated  Legendre  polynomials  of  both  kinds 
and  their  derivatives  are  computed  [6  J  and  stored  in  the 
machine.   Then  the  Fourier  coefficients  are  found  directly 
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from  eqo  (5»21),  and  the  systems  of  equations  given  by 
is, 22)   must  be  solvedo   Elimination  methods  do  not  take 
advantage  of  the  great  number  of  zeros  in  the  matrix  of 
the  coefficients  and  so  are  very  time  consuming.   They 
also  give  rise  to  very  large  errors,  in  some  cases,  due 
to  loss  of  significance,   A  method  suggested  by 
S.  Schechter  [?]  is  used  and  is  found  to  be  very  accurate. 

(    *     \  /    \ 

The  computation  of   5W     and   5W      is  entirely 

(e) 
straightforward.   For   5W      it  was  observed  that  the 

computation  of  the  double  sum  given  in  eq,  (5o26)  was 

the  most  time  consuming  part  of  the  whole  code^   Instead 

of  simming  from   -N   to   N-1,   we  tried  summing  from 

-N«   to   -N»-l   where   N«  <  No   The  omitted  terms  were 

found  to  make  a  very  small  contribution  and  the  code 

could  be  speeded  up  in  this  way  by  a  factor  of  as  much 

as  two.   Our  results  show  that  fairly  good  Indications 

of  stability  can  be  obtained  with  N  =  8   and  K'  =  6. 

However,  in  regions  where   ^    is  small,  the  matrix 

elements  fall  off  very  slowly  with  increasing  n  and  y\L 

as  can  be  seen  by  an  examination  of  the  equations,  and 

the  matrix  segment  with   N  =  N'  =19,   the  m^axlmum 

permitted  by  the  code,  may  not  be  large  enough.   This  is 

not  a  very  serious  limitation  since  these  are  not  regions 

of  physical  Interest.   Once  the  matrix  of      5W   is  obtained 

the  eigenvalues  are  computed  [8 J  by  rotating  the  matrix 

to  tridiagonal  form  and  then  using  a  Sturm  sequence  of 
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polynomials  to  isolate  the  roots » 

The  computations  show  that  for  a  given  configuration 
there  is  a  value  of   N   (and  N* )   such  that  further 
increase  in  size  does  not  change  the  eigenvalues  — 
the  added  terms  are  so  small  that  the  added  eigenvalues 
are  computed  to  be  zeroc   If  N   (and  N' )   are  taken  to 
be  less  than  optimum,  it  has  been  observed  that  the 
resulting  eigenvalues  vxill  be  changed  in  the  direction 
to  make  the  system  appear  less  stable.   Thus  the  machine 
computations  give  results  which  are  only  a 
lower  bound  on  the  complete  region  of  stability. 

Per  a  given  m,   it  can  be  shown  that   5W   depends 
on  four  parameters,  the  aspect  ratio   A ;   the  ratio 
between  the  conductor  radius  and  the  plasma  radius  denoted 

by   A.  s   9-nd  two  further  parameters  which  can  be  chosen 

2 

as   B.J   the  ratio  of  the  '''     \-^e"'  magnetic  pressure  of 

the  toroidal  field  inside  t    _  ^asma  to  the  gas  pressure, 
and  B^,   the  ratio  of  the  "a\'erage"  magnetic  pressure 

of  the  toroidal  .  '     in  the  vacuijm  to  the  gas  pressure. 

„2  Ap 

Bg      and     B^'      are    defined  by 

B~   =  b^/(2a2p)  sf  =  b2/(2a2p}   .  (6.1) 

Thus    in   computing      6W     we    can  normalize   both      a      and      p 
to   unity. 

Now    ir:  r   to   make    comparison  with  the    cylinder  we 
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define  similarly  dimensionless  field  quantities  for 
the  cylinder 

b^=b^/2p  b^=b^/2p  .  (6o2) 

It    is    convenient   to  use   the    same   units   as   does   Tayler. 
He    defines 


b.    -b^/B^  b^  =  b^/B^  {6.3) 


and   substitution  of    (6o3)    into    (Lt.,2)    gives 

1   +  b^    -   b^     =     2p/Bw     .  (6„l^.) 

Since  the  right  hand  side  is  positive^  we   have  the 
restriction 

1  +  b^  >  b?  o  (6„5) 

e  —   1 

A 

To  relate   b.   and  b    to   b    and  b^  we  set 
1        e       c        c 

(6o2)  into  (I4.02)  giving  an  expression  for  ^^Z  V     which 
is  substituted  into  (6oI).)  .   Then  after  some  reduction 
we  get 

A 


b^  =  b^  yi/{l+b^-b^)     b^  =  b^  Vi/(l+b^-b^)»  (6c6) 

From  (2o8)  we  have 

1  +  G^  (b!-B^)   =  bI/  2p  (6c7) 

where  the  geometric  factor   G  =  (cosh  K^   -  cos  0)/sinh  )^o 
In  the  limit  as   >7   —>  00  »   G  — 5>  1  but  even  for  moderate 
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values    of     />?    ,       G      Is    close    to   unityo      The   right    hand   side 
of    (60 7)  must   be   positive    so  we   have    approximately 

1   +  B?     >      B^  „  (6.8) 

1     —        e 

For  comparisons  between  cylinder  and  torus,  we  choose 

critical  values  between  stability  and  instability  of 

b    and  b.   and  use  (6,6)  to  get  the  corresponding 

values  of  b    and  b  .   These  are  then  Identified 
c        c 

?        2 
with  Bo   and  B  <, 
1        e 

An  examination  of  equations  for   5W  (.S'^^-SoZJ ) 
would  lead  one  to  expect  that  configurations  which  are 

stable  for  ra  =•'  1  vfill  be  stable  for  higher  values  of 

(i^  (e ) 

mo   5W   ^   and   5W'  "   are  derived  from  positive 

definite  expressions  and  so  must  have  eigenvalues  which 

are  all  posit ive.   These  terms  vary  as  m   »   The  surface 

term  (3, 2?)  gives  a  negative  contribution  to   5W  but 

V8J»ies   as  m  ^o   Eigenvalues  of  the   5¥  matrix  have 

been  calculated  for  m  =  2   and  m  =  3o   The  results 

show  that  if  m  =  1   is  stable  m  -"-  3   is  always  stable* 

In  a  few  cases  m  =  2  is  not  quite  stable^  but  so  close 

to  stability  that  only  a  slight  shift  in  one  of  the 

parameters  would  be  sufficient  to  achieve  stability* 

In  general  it  can  be  stated  that  if  m  =  1   is  stable 

all  higher  modes  are  also  stable.   The  case  m  =  0   has 

been  discussed  above ^   The  results  quoted  below  are  for 

m  =  1   onlyc. 
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Regions  of  stability  have  been  mapped  out  by  choosing 
B.  ,  B    (satisfying  6.8)  and  yl  =  R.  /  R    and  then 

X     '         e  d         a  J.       I'o 

determining  what  value  of  /\     will  give  a  stable 
configuration.   It  has  been  found  that  a  stable 
configuration  can  be  achieved  for  very  large  ranges 
of  the  parameters.   Roughly,  for  -A   as  large  as 
5   a  stable  configuration  exists  for  all   B.   and   B   <  10. 
If   B.   and   B    do  not  differ  by  more  than  ten  or  twenty 
percent  the  ranges  of  these  parameters  for  which  stability 
exists  is  extended  up  to  100  or  more.   For  a  given  triplet 
B.  ,  B    and  j\      in  these  regions  there  is  a  critical 

J-      o 

A.         such   that    for    all      /\    >      ^        the    conl'igur ation   is 
—  c  c 

unstable   while    for    all      ^     <      A         the    ccnf iguiration   is 

stable.    ^      It    should  be    noted   that    In   general      Z\         is 

'='  c 

small,  i.e.  lo3  or  less.   To  achieve  stability  I'or  an 

aspect  ratio  as  large  as  3*0  requires  that   B.   ti;'  B 

(  s ) 
(to  make  the  negative  contribution  I'rom   5W     very 

small)  and  that  J'y      Pe  close  to  unity. 

The  results  can  be  presented  in  another  way  which 

gives  some  comparisons  between  torus  and  cylinder.   From 

eq.  (I4-.I5)  values  of   b   ,  b    ana  /\.    are  determined 

such  that   5W  =  0   i.ec  the  cylinder  is  in  neutral 

equilibrium  for  ,gne  value  ol'   k  but  stable  for  all 

other   kc   This  value  of  k  (=k  )  determines  a  wave  length. 

If  a  cylinaer  of  just  this  length  is  bent  into  a  torus 
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then  the  stability  or  instability  of  this  toroidal 

configuration  can  be  taken  as  an  indication  of  the 

effect  of  the  curvature »   In  making  the  correspondence 

between  cylinder  and  torus,  the  magnetic  fields  are 

identified  as  explained  above.   The  aspect  ratio 

A   =  cosh  >!-,   is  computed  by  identifying   ^^n^i^ 

with   (sinh  1-,  )~    «   The  computations  show  that  in 

most  cases  the  toroidal  configurations  obtained  in 

this  way  are  stable,  (Exceptions  lie  in  the  regions 

where   B.   and   B    are  large  and  not  nearly  the  same.) 
1         e 

It  should  be  noted,  however,  that  the  aspect  ratios 
in  these  cases  are  all  small  --  less  than  1«25  for 
yi  =  3oO   and  less  than  lo75   for   A  "  l»5o 

These  results  indicate  that  the  effect  of  curvature 
is  stabilizing  when  the  torus  is  related  to  the  cylinder 
in  the  manner  discussed  in  the  preceding  paragraph. 
On  the  other  hand,  a  thin  torus,  that  is<,  one  with  a 
large  aspect  ratio^  is  never  stable  and  in  fact  this 
study  indicates  that  it  is  not  possible  to  obtain  a 
stable  plasma  in  a  toroidal  gas  discharge  which  has 
both  a  reasonable  aspect  ratio   (A  >^  3 )   and  a  value 
of  J\     large  enough  to  compress  the  gas  well  away 
from  the  conducting  walls » 
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